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Abstract. We solve the two-particle s-wave scattering problem for ultracold atom 
gases confined in arbitrary quasi-one-dimensional trapping potentials, allowing for two 
different atom species. As a consequence, the center-of-mass and relative degrees of 
freedom do not factorize. We derive bound-state solutions and obtain the general 
scattering solution, which exhibits several resonances in the ID scattering length 
induced by the confinement. We apply our formalism to two experimentally relevant 
cases: (i) interspecies scattering in a two-species mixture, and (ii) the two-body 
problem for a single species in a non-parabolic trap. 
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1. Introduction 

A strongly interacting ultracold atom gas displays interesting features of a correlated 
quantum many-body system when its dynamics is confined to one dimension pQ. The 
presence of a transverse confining potential has been shown to induce characteristic 
resonances in the coupling constant of the two-particle s-wave scattering process 
which have become known as confinement- induced resonances (CIR). The existence of 
the CIR has been revealed under the simplifying assumption of a transverse parabolic 
confinement potential with length scale a± and for the case that the two scattering 
atoms belong to the same species El Hj- In this case, the center-of-mass (COM) 
and relative coordinates of the two particles can be separated, allowing to factorize the 
problem into single-particle problems. At low temperatures, only the COM ground state 
is occupied, the decoupled COM motion can be disregarded, and the two-body problem 
can be solved exactly within the pseudopotential approximation. The result is that 
there is exactly one bound state for any 3D scattering length a. In the limit of small 
binding energy, the particles are tightly bound in the lowest-energy transverse state and 
form a very elongated dimer. The appearance of such a bound state is purely due to 
the confinement, since for a < no dimer is formed in free space. In the opposite limit 
of large binding energies, the dimer becomes spherically symmetric. In this regime, the 
confinement is not effective, and the free-space result is recovered. Moreover, a unitary 
equivalence exists between the Hamiltonian and its projection onto those channels which 
are perpendicular to the one with lowest energy. As a consequence, to each bound state 
corresponds a bound state of the closed channels, which then causes the CIR |3J. It 
occurs at a universal value of the ratio aj_/a = C = —£(1/2, 1) ~ 1.4603, where ( is the 
Hurvitz zeta function. The influence of the CIR has also been studied for the three- 
body [3] and the four-body problem jE] in the presence of confinement. In particular, 
the solution of the four-body problem completely determines the corresponding quasi- 
lD many-body BCS-BEC crossover phenomenon Recently, the existence of the 
confinement-induced molecular bound state in a quasi-lD Fermionic 40 K atom gas 
confined in an optical trap has been reported [7] . By using rf spectroscopy, the binding 
energy of the dimer has been measured as a function of the scattering length, with 
quantitative agreement to the results of Ref. [HJ. However, the existence of the CIR in 
the scattering states remains to be observed. 

Although the analytical results for the parabolic confinement are instructive, 
realistic traps for matter waves frequently have non-linear potential forms, see for 
instance Ref. [Hj for a particular example of a trap on the nanometer scale. To give 
another example, for the problem of tunneling of a macroscopic number of ultracold 
atoms between two stable states of a trapping potential, the nonlinearity clearly is 
crucial. Hence generalization to the non-parabolic case is desirable and provided in this 
work. In addition, we consider traps with two different species of atoms. Note that 
sympathetic cooling techniques require to study this case. Different trap frequencies 
may arise for different atom species, e.g., because of different atom masses or different 
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magnetic quantum numbers. Here we obtain general expressions for the bound-state 
energies and scattering resonances when the COM and the relative degrees of freedom 
do not decouple anymore. In the parabolic limit and for intraspecies scattering, we 
recover well-known results E]. For the general case, we show that more than one 
CIR may appear, and that it depends on the symmetry properties of the confining 
potential how many resonances occur. We apply our formalism to two experimentally 
relevant cases: (i) interspecies scattering in a two-species mixture of quantum degenerate 
Bose and Fermi gases in an optical trap, and (ii) a single species cloud in a magnetic 
trap, taking into account non-parabolic corrections due to a longitudinal magnetic field 
suppressing Majorana spin flips. 

As we will discuss below in more detail, the CIR has a close similarity to the well- 
known Feshbach resonance which arises if the Hilbert space can be divided into 
open and closed channels coupled together by a short-range interaction. Due to this 
small but finite coupling, two incoming particles initially in the open channel visit the 
closed channels during the scattering process. If a bound state with energy close to the 
continuum threshold exists, such a process is highly enhanced and a resonance results. 

The paper is organized as follows: In Sec. we present the general formalism. 
Section E] presents the bound-state solution, while Sec. 0] contains the analysis for the 
scattering solutions, including the analogy to Feshbach resonances. In Sec. El we discuss 
the special case of harmonic confinement, and in Sec. El a particular example of a non- 
parabolic confinement is illustrated. Finally, we conclude in Sec. Technical details 
have been delegated to two Appendices. We set h = 1 throughout this paper. 

2. The two-body problem 

Let us consider the general case of two different atomic species with mass mi and m%. We 
denote the particle coordinates by x, = (xj^, Zj) and their momenta by pj = (pj^piu). 
Different atoms may experience a different transversal confinement potential Vi(x.± t i). 
For ultracold atoms, only low-energy s-wave scattering is relevant, and the interaction 
between unlike atoms (and similarly, also the interaction between the same atoms) 
can be described by a Fermi-Huang pseudopotential V(|xi — x 2 |). Then the relevant 
Hamiltonian for two different atoms is given by 



where /i = 777.1777.2/ (mi + 7772) is the reduced mass and a the 3D scattering length. This 
allows to characterize the two-body interaction by the parameter a only. The validity 
of the pseudopotential approach has been verified numerically for finite-range potentials 
in Ref. For further convenience, we transform to the relative/COM coordinates 




(1) 



The pseudopotential has the standard form 




(2) 
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and momenta given by r = (r±,z), R = (Rj_, Z) and p = (p±,p\\), P = (Pj_,Pn), 
respectively. This can be done by the canonical transformation 

/R\ ( m l m 2 \ / Xi\ 

1 M -M 

M M M 

\ m% —mi ) 

where M = m\ + m-i. Since the confinement is assumed to be purely transversal, 
the longitudinal COM coordinate Z is free and decouples from the other degrees of 
freedom. Hence we eliminate it by transforming into the longitudinal COM rest frame, 
where the state of the system is determined by the set of coordinates (xj^i, Xj_ j2 , z) 
or, alternatively, by (Rj_, r) = (Rj_, rj_, z). The transformed Hamiltonian takes the form 

H = H\\+Hj_, 1 + H_ L ,2 + V, (4) 

where 



m = P X 



pi, 



(5) 



2\x ~ ' 2mi 

For a more compact notation, we introduce the non-interacting Hamiltonian H = H—V 
and denote its eigenstates by 



|fc,A 1 ,A 2 ) = e- ifc2 ^ 1 i ) (x ± , 1 )^(x ±i2 ), 

where ip® are single-particle eigenstates of H ±i for eigenvalue E^}. 
the two-particle Schrodinger equation is given by 

(H - E) *(R ± , r) = -y(r)¥(Rj_, r) . (7) 

The pseudopotential (J2J) can be enforced by the Bethe-Peierls boundary condition 

/(R±) /- r 



,( 2 )/ 



(6) 

Corresp ondingly, 



* R 



0) 



Airr v a, 
leading to the inhomogeneous Schrodinger equation 

/(R±) 



(H -E)V(R ± ,r) 



2(JL 



-8(r). 



(8) 



(9) 



The solution of this equation can be formally obtained in terms of a solution of 
the homogeneous Schrodinger equation, (Hq — E)^q = 0, and the Green's function 
Ge = {Hq — E)~ , 



*(Rj 



^o(R±,r) + J dR^G^(R ± ,r;R^0)- 



2fi 



(10) 



To determine /(Rj_), we substitute Eq. (jTUj) into Eq. (jHJ) and find the integral equation 

/(R±) _ , T , fT) n \ , f JX3 / /• /T3 D' \f/T)' ^ ^-q 



4-7ra 



^o(R±,0) + y dRlCs(R±,Rl)/(R^), 



where we have defined the regularized integral kernel 

Ce(R-U R-l) = lim ^ (G S (R ± , r; R' ± , 0) - 8(R X - R' ± 



27rr/ 



(12) 
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In Eqs. (JTUJ) and \& can be expressed as a superposition of single-particle 

eigenstates \k, Ai,A 2 ) with 

| + B <.) + s p) = B . (13) 

We refer to the set of states with the same transverse occupation numbers Aj but 
arbitrary longitudinal relative momentum as a scattering channel or, simply, channel. 
Each channel has a minimum energy given by + E x . Since the interaction is 
short-ranged, only states fulfilling Eq. (fTHJ) appear in the asymptotic solution. For 
each open channel, E > E^' + E^\ such that there are (at least) two such states 
having opposite momenta. For E just above Eq + E^\ there exists one open channel 
only. The corresponding solution given by Eq. ()1()J) describes the scattering of two 
particles initially occupying the transverse ground-state. During the scattering process, 
the particles populate closed channels, but afterwards return into the single available 
open channel (quasi-lD picture). For E < E^ + E^\ all channels are closed and only 
bound-state solutions are possible. These are given by Eq. (JTUJl with \I/ (R_|_, r) = 0. In 
the following, we consider both classes of solutions in more detail. 



3. Bound-state solutions 



Let us consider the situation when all channels are closed and only bound states may 
occur. We define the binding energy of the bound states as 

E B = E - E > , (14) 

where E = E^ + E^ is the ground-state energy of H . To find bound states, we 
diagonalize the operator £ E (R_i_, Rj_) defined in Eq. (fT2j) . where Eq. (fTTj) yields the 
condition 

- = J dR' ± C E (R±, R/l)/(R1) • (15) 

For given a, bound states with binding energy Eb = Eq — E follow as solution 
of this eigenvalue problem. The bound-state wave function follows by inserting the 
corresponding eigenvector f(R±) into Eq. (fTUj) with \E f (R_ L ,r) = 0. In order to find 
a representation of Cs(R±,Rj_) allowing for straightforward analytical or numerical 
diagonalization, we use 

POO 

G E (R ± , r; R' ± , 0) = / dte Et G t (R L ,r;R' ± ,0) , (16) 
Jo 

with the imaginary-time evolution operator 

G t (R ± ,r;R' ± ,0) = (Rj_, r| exp[-#o*]|R±, 0) (17) 

for H Q . The time evolution operator exp[—H t] can be factorized into the product 
exp[— H\\t] expj— Hi t ±t] exp[— H2,sj\- The corresponding factors in G t are 

(,|exp[-^t]|/) = (^) 1/2 e-(-')V- (18) 
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for the relative longitudinal coordinates and 

(x^lexphi^lx^) = Y^ e ~ E ^x(^x\<±) (19) 

A 

for the transverse coordinates (the bar denotes complex conjugation). Thus Gt can be 
expressed in terms of the set of coordinates (x_l 1 ,xj_ 2 , z) as 



G? t (R ± ,r;Rl,0) = J^'^ II E e ~ 4 M° (^MVx,) • (20) 

' 71 i=l,2 A 

This equation illustrates that for large imaginary times the integrand in Eq. ([THJ) decays 
as exp[— Est}. Notice that this representation is valid for Eb > 0. By using 



oc 



dtl-?-) e~ r ^ 2t (21) 
2vrr L \2ntJ v ' 



we find 



r°° dt r ^ / ii \ 3/ 2 i 

Ce{R ± ,B! ± ) = I ^[e^(R ± ,0;Rl,0)-(^j 5(R ± -R' ± )] .(22) 

To show that the integral in Eq. (|22j) converges also for small t, we expand 
Gi(R_i_, 0; Rj_, 0) with respect to t, see Appendix A. We find 



hmG 4 (R ± ,0;R^,0) = S(R ± - R' ± ) - ( ^ 



2nJ 



x 



P 



2 



+ ^i(R±) + ^ 2 (Rj 



(23) 



_2M 

Thus (e can be regarded as a regular operator acting on the space C 2 of square-integrable 
functions. We note in passing that if the two single-particle transverse Hamiltonians 
H±j commute with the angular momentum operators L z , then also (e commutes with 
L z . This follows by observing that in this case we can choose for the eigenbasis {4\ } a 
set of eigenvectors of L z , and the product of two eigenvectors of L z is still an eigenvector 
of L z . Hence Ce(R_u Rj_) can t> e written as a sum of projectors onto states with definite 
angular momentum. A similar conclusion can be drawn regarding parity symmetry, 
when considering non-cylindrical confining potentials that obey this symmetry. 
Using the transverse non-interacting ground state, 

V*(R±, O = 4 1] (r± + —r±) 4 2) (*± - —t±) , (24) 



mi J \ m 2 

provided the overlap integral f c/R / J _^o(R / J _, 0)/(R' ± ) ^ 0, the integrand in Eq. (|TU|) 
decays as exp[— Est]. This defines a spatial scale a# for the longitudinal size of the 
corresponding bound state, as = 1/^/JlE~b. For large Eb, clb is small and we have 
very tight pairs. This constitutes the dimer limit. On the other hand, for small Eb, 
atom pairs are very elongated. This regime is termed BCS limit. In the following, we 
investigate both limits in greater detail. 
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3.1. Dimer limit 

For large binding energies, the atom- atom interaction dominates over the confinement. 
Due to the exponential factors in Eq. f!22L only small imaginary times contribute 
significantly to the integral, and we can substitute G t with the short-time expansion 
(J23J) as derived in the Appendix A, yielding 

-e Et t- 1/2 (R ± \H ± \R' ± }) . (25) 

Hence, the operator ( B now shares eigenfunctions with H± = Pj_/2M+Vi(R_l)+V2(Rj_). 
For (identical) parabolic confinement potentials, H±_ is exactly the decoupled COM 
Hamiltonian. Let us denote the eigenfunctions and eigenenergies of H± as 0a(R±) and 
E^\ respectively. Substituting 4>\(R±) into Eq. ()15|) yields after some algebra 

+ (26) 



Ana 47r \ 2\E\ I 4n 

In the second relation, we have used Eq. (jl4j) . From this, we directly obtain the binding 
energy in the dimer limit a — > + as 

1 



(27) 



2/xa 2 ' 

which coincides with the result obtained in free (3D) space without confinement. 
3.2. BCS limit 

The scattering channel with lowest energy, corresponding to the transverse non- 
interacting ground state ip , opens at the energy threshold E = E . For E B — > + , 
as the energy approaches this threshold, the term with Ai = A2 = dominates in Eq. 
(I2UJ) . and yields in Eq. the contribution 



— ^o(R±,0$ o (R±,0), (28) 
8fiE B 

which diverges for E B — > + . All other channels are still closed at E = Eq and give 
finite contributions in Eq. (|20j) . This observation suggests a useful separation of the 
total Hilbert space into a part 7i corresponding to the open channel (or lowest-energy 
scattering channel) and a part 7i e perpendicular to that. With this separation, terms 
yielding a finite contribution at E B — > + can be summarized in the Green's function 

G t (R ± ,r;R' ± ,0) = (R ± , r| exp[-H t]\R' ± , 0) , (29) 

where Hq is the projection of Hq onto the Hilbert subspace H, e . We then define a new 
integral kernel, 

r°° dt r / u \ 3/2 1 

Cb(R±,R1) = J -[e^(R ± ,0;R' ± ,0)-(^j 5(R ± -R' ± )] , (30) 

which is also well-defined for energies above the threshold E = E . 
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For small Eb, Eq. (J15)) is most conveniently solved by expanding /(Rj_) in an 
orthonormal basis \j ) according to 



/>=E£U>> fs= fdR±(j\R±)f(B.±), 

3 J 



(31) 



where the basis state | ) corresponds to 

(Rl|0> = c^ (R±,0), (32) 

with normalization constant c. Although ifjo(FL±, r) is a normalized element of the two- 
particle Hilbert space, this does not imply that i/jq(R,j_,0) is an element of the COM 
Hilbert space with norm unity. In fact, the normalization constant c has to be computed 
explicitly and generally depends on the particular confinement. In this basis, Eq. (|15j) 
assumes the compact form 



|0) is an approximate eigenstate for small Eb, since all the matrix elements are finite 
apart from (0|£e|0) which diverges according to 



< 0|&|0> ^/8^? + < 0&|0> - (34) 



Substituting this in Eq. (|3lj|) yields 



' ' ' 1 + (0\Ce\0) ■ (35) 



Ana V 8^iEb c 2 

Neglecting the last term, the relation for the binding energy Eb is solved in the BCS 
limit a — » 0~, 

2a 2 n 2 , . 

E B w — — • 36 

/iC 4 

4. Scattering solutions 

In this section, we focus on scattering solutions at low energies E slightly above Eq, 
where exactly one transverse channel is open. Then the incoming state is given by 

M/o = e lfc2 4 1) (x li± )4 2) (x 2 , ± ), (37) 

which describes two incoming atoms with (small) relative longitudinal momentum 
k = ^JlmiE — E ) in the (transverse) single-particle ground states ip^ and ^\ 
respectively. 
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4-1. One- dimensional scattering length a^o 

As done in Sec. El we split off the contribution from the open channel, 

dte Et G t (R ± ,r;-R' x ,0), (38) 

where G t (Rj_, z; Rj_, 0) is the Green's function restricted to TC e , which is well-defined also 
above E . Inserting Eq. (J3~8~j) into Eq. (fTUj) yields for \z\ — > 00 the standard scattering 
solution, 

tf(R,r) = ^o(R±,rx)(e^ + / e (fc)e ifc l 2 l), (39) 
with scattering amplitude 

fe(k) = ^J dR' ± MR' ± , 0)/(R' ± ) , (40) 

whereas for short distances, also the term J dR! ± J Q °° dt e Et G t (R±, r; Rj_, 0)/(R'j_) 
appears in the scattering solution. Since the energy is well below the continuum 
threshold for the closed channels, this must be regarded as a sum over localized states. 
Enforcing the boundary condition (|HJ) then leads to an integral equation for f(R±), 

~^T = / rfR/ ^(R±,Rl)/(RY) 

+ ^ (R±, 0) + #o( ^' 0) J dK' ± MK' ± , o)/(Rl) . (41) 

This integral equation is most conveniently solved by again expanding /(Rj_) in the 
orthonormal basis {\j)} introduced in the previous section. Thereby, we can express 
Eq. ()41|) in compact notation, 

which is formally solved by 

1/)= , ~)!, C J Ce + ^-) |0>. (43) 



1 — if (kaii)) V ^ 7ra 
The parameter a^ follows in the form 

2c 2 

a iD = = . (44) 

(0|[C s + l/(4 7 ra)]- 1 |0) 

From Eq. ()40|) . f e {k) = —1/(1 + ikany), which allows to identify aiD with the ID 
scattering length. Having introduced this parameter, the ID atom-atom interaction 
potential can then be written in an effective form according to 

V m (z,z')=gi D 8(z-z'), (45) 

with interaction strength giD = —l/(fian)) [2 . For very low energies, k — > 0, we can now 
formally set E = Eq in Eq. (}4*4*j) . For a confining trap, Qe q is an Hermitian operator with 
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discrete spectrum {A n } and eigenvectors |e n ), which eventually have to be determined 
for the particular Hamiltonian. Thus we find 

a ~ 1 V K°l^>l 2 (m) 
^-^^ An + i /(4vra) - ^ 

This result has interesting consequences for the two-body interaction. The denominator 
can become singular for particular values of a, thereby generating a CIR. Every 
eigenvalue \ n corresponds to a different CIR, unless the overlap (0|e n ) vanishes due 
to some underlying symmetry of the Hamiltonian. We anticipate that for identical 
parabolic confinement potentials, the decoupling of the COM motion implies that only 
one resonance is permitted. For confining potentials with cylindrical symmetry, there 
is a resonance for each eigenvector of (e with zero angular momentum. For confining 
potentials obeying parity symmetry, the eigenstates |e n ) must be even. These two 
symmetries allow in principle for infinitely many resonances. In practice, however, only 
few of them can be resolved because the resonances become increasingly sharper when 
|(0|e ra )| 2 — > 0, making them difficult to detect. 

4-2. Interpretation of the CIR as Feshbach resonances 

A very simple and illuminating analysis, similar to that for standard Feshbach resonances 
jHj, is also possible for the CIR. The two-particle Schrodinger equation can be rewritten 
as an effective Schrodinger equation for the scattering states in the open channel, 
(E — H c ff)V\ty) = 0, with the effective Hamiltonian 

H eS = #opcn + VHM MHV ' ( 47 ) 

J-J closed 

Here, H opcn = VHV and H dosed = A4HA4, where V and M. are projectors to open and 
closed channels, respectively. This equation can be expressed in terms of the closed- 
channel eigenstates 

H eS = H open + VHY\ } (48) 

with i/ c i ose d 1 3> n ) = E n \<& n ). This implies that a Feshbach-like resonance is possible 
at zero momentum if two conditions are fulfilled. First, there exists a solution of 
(i?o — inclosed) | = 0, i.e., |<&) is a bound state of inclosed with energy E = E . Second, 
|$) must be coupled to the open channel, VH\Q) ^ 0. 

Within the pseudopotential approximation, the equation (Eq — -ffcioscd)!^) = is 
solved in terms of the Green's function 

POO 

MG Eo (R ± ,t;R' ± ,0)M= dte Eot G t (R ± ,r;R' ± ,0) (49) 

Jo 

by the state 

$(Ri,r) = J dR' L MG Eo (RL,r;R' L ,0)M^p- , (50) 
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together with boundary condition 



Aixr \ a 
This leads to the eigenvalue equation 

-£ = &„!/), (52) 

which is solved by the eigenvectors \e n ) introduced above. This yields a = — l/(47rA n ), 
implying that there is a bound state |$) of H dosed with energy equal to the energy of 
the incoming wave, corresponding to the resonances found in the previous subsection. 
The CIR is then in complete analogy to a zero-momentum Feshbach resonance. Due to 
the small but finite coupling to the closed channels, two incoming particles initially in 
the open channel visit the closed channels during the scattering process. This process 
is strongly intensified when a bound-state exists whose energy is close to the continuum 
threshold. Then, a scattering resonance results. Note that such a bound state can be 
occupied only virtually by two particles during the scattering process. Hence from now 
on we will refer to such a bound state as a virtual bound state. 

It is also possible to recover the overlap condition (0|e n ) ^ in this framework. In 

fact, 

PH$(R ± , r) = PV(t)$(R ± , r) = -V^^-5(r) = -^o(R±, t ± )S(z)^- , (53) 

since |$) fulfills Eq. (15 1|) with /(R) = (Rj_|e n ). Hence, the two overlap conditions 

VH\$) ^ <£> (0|e n ) ^ (54) 

are equivalent. When they are not fulfilled, there exists a virtual bound state with 
energy E , but it is not coupled to the incoming wave. 

5. Special case of harmonic confinement 

In the previous sections, we have formulated the theory for a general confining potential 
and for two different atomic species. As a simple illustration, we now consider the case 
of harmonic confinement, V^(xj) = m i uj 2 Ti 2 _ L /2. In COM and relative coordinates, 

1 1 fa 2 u 2 \ 

Konf(R±, r±) = - (m^f + m 2 u)l) |RJ 2 + - — ui\ + — \r ± \ 2 

+ /i (cof - uol) r± R± . (55) 

In general, the COM and the relative coordinate do not decouple, and in order to find 
the scattering and bound-state solutions, we have to follow the procedure outlined in the 
previous sections. To that end, we label the single-particle transverse states by quantum 
numbers A = {m, n}, where m is the integer angular momentum and n the integer radial 
quantum number. The eigenenergies and -states of the 2D harmonic oscillator 

F {i) -nc ih {i) - 1 ih f X± 
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{rriiUJi) 1/2 ,i 



1,2, can be expressed in terms of the 



with the oscillator lengths cij 
quantities 

e n , m = 2n + \m\ + 1 and ip n ,m(^±) = e im ^Rn, m (|xjJ), 

where 



R n ,m{p) 



nl 



1/2 



^7T \(n+ |^|)!, 

with L„ (a;) being the standard Laguerre polynomials. A convenient choice for the 
orthonormal basis \j) introduced in Eq. (|ST|) is then given by 



(Rjj) = (Rjm,n) = — W 

fl M 



R, 



(56) 



with the length scale au = (tojCJi + m 2 w 2 ) 1//2 . In particular, we find for |0) = |0,0) 
that (Rj.|0) fulfills Eq. (J22D with c = v /ia 1 a 2 /a M . 

The single-particle imaginary-time propagator for a 2D harmonic oscillator with 
length scale a$ and frequency u is given by 



vra 2 , 1 - e" 2aJ * 



exp 



xj + x'j 
2a 2 



coth(u;£) 



a 2 , sinh (u;i) 
(57) 



Inserting this into Eq. (|20j) with x_l, 
G t (R ± ,0;Ri,0) = 



R 



R'i and z = 0, we find 



/x 0(1-/3) 



-UJ\t 



-L02t 



2jit 7r 2 a^ 1 - e" 2 "** 1 



R 2 +R 



/2 



2rt 2 



/(*) 



g— 2i<J2t 

R i ■ R' 



x exp 

where we have introduced (3 = a 2 M /a\ and 

f{t) = /3coth(wit) + (1 - /3) coth(w 2 *) , 
= /5sinh _1 (u;it) + (1 - 0) sinh -1 ^) . 



a 



g(t) 



M 



(5? 



(59) 



In order to compute explicitly the operators (e and we still have to project onto the 
discrete basis {\j)} and to perform the imaginary-time integral for each matrix element. 
In general, this cannot be achieved analytically, and one has to resort to a numerical 
evaluation. Only for u)\ = u 2 , a complete analytical solution is possible. Since the COM 
degrees of freedom separate, this solution is a trivial extension of Ref. j2]. Nonetheless, 
along with the general analysis of the previous section, it provides a physical picture for 
weak interaction between the COM and the relative degrees of freedom. 



5. 1 . Identical frequencies 

For lo\ = 0J2 = uj, the COM and relative coordinates separate, H = H Te \ + Hqom, with 



"■' = I + 2^ + V{T) 



COM= 2M + 2 M " 2Ri 



(60) 
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In this case, we can consider the two-particle system being (asymptotically) in the 
ground state of the decoupled COM Hamiltonian, and just solve the relative problem 
[21 IS] - Moreover, with f(t) = coth(u;£) and g(t) = sinh _1 (u;t), the Green's function (|58j) 
simplifies to 



G t (Rx,0;Ri,0) = J\ ^L ^^^m (^) i>n, m {^A -(61) 

V 2nt 7ra 2 M 1-e 2uJt </ \a M J \a M J 

In this case, \n,m) is an eigenstate of the decoupled Hamiltonian H C om, and describes 
the COM motion also for finite r. Moreover, au = (Mo;)" 1 / 2 and a„ = a M /{p(\-P)) = 
(fiuj)' 1 / 2 are the characteristic lengths associated with Hqom and if re i, respectively. 
Inserting Eq. (j61|) into Eq. (|22j) and rescaling t by 2u, we obtain 

s ^\n,m)(n,m\ f°° dt / e -^ME)t x \ 
^ " JT n / TZ7\T/2 1 _ p -t ~ 7 ' (62) 



4vra M 7 (7ri)V2 ^ i _ e -t t 

with Q ntm (E) = (1 + e niTO — E/u)/2. The integral on the rhs of Eq. (jn2*|) is related to 
the integral representation of the Hurvitz zeta function ^(1/2, 0„ iTO ) |51 ITU]. 

5.1.1. Bound states The condition given in Eq. ()15|) for a bound state with transverse 
configuration \n, m) translates into 

1 ~ \ a. 



C(^n, m ) =-^f. (63) 

The zeta function is monotonic, and has the asymptotic scaling behavior 

CQ, n<l) «JT 1/2 , cQ,^»l)^-2v / fi. (64) 

Inverting Eq. ([63)1 . we recover the bound-state energy found in Ref. The 
corresponding result is plotted in Fig.[T] As an immediate consequence of the decoupling 
of the COM degrees of freedom, the eA-fold degenerate energies corresponding to excited 
transverse configurations follow from the COM transverse ground state by a shift along 
the ordinate in steps of u>. This is indicated by the dotted curves in Fig. [TJ Notice 
that for energies above E = 2u>, corresponding to E B = 2uQ 00 (E) < 0, there exists 
an open channel, but the solutions associated with COM excited states are orthogonal 
to it. For this reason, the relevant condition for a bound state to exist with transverse 
configuration \n, m) is fl n m (E) > 0. /,From the scaling behaviors in Eq. ([64)1 . we find 
the limiting behaviors of the energy of the bound state at \a^/a\ <C 1 as 

E B ,n,m ~ 77— T for « > , 

2o? 

E B ,n,m ~ — 7 for a < , (65) 

see Eqs. (|2Tj) and (|3^j) . with c = y/na^ and E Bn ^ m = u>Q n>m . Hence, in this highly 
degenerate case, there is exactly one bound state for each transverse configuration and 
each scattering length a. 
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5.1.2. Scattering states In order to identify resonant bound states of the closed channel, 
and the corresponding zero-momentum CIR, we subtract the contribution of the lowest- 
energy scattering channel in Eq. (JH2J), and obtain 



with Q fi(E) = Qq^(E) + 1 and VL nm (E) = Q nm (E) for n + |m| > 0. Hence the 
curve corresponding to the COM ground state is shifted vertically by 2uj and coincides 
with the curve corresponding to the excited states |1,0), |0,2) and |0, — 2). Moreover, 
the coupling condition in Eq. (J53J) becomes (0,0|n, m) ^ 0, and is fulfilled only for 
n = m = 0. Though there are in principle infinitely many closed-channel bound states 
with energy 2u (one for each curve), only one scattering resonance exists, since only one 
of them is coupled to the incoming scattering wave. Inserting Eq. (|66|) into Eq. (|4*Bj) we 
recover for the ID interaction strength giD the well known result j2] 



5.1.3. Physical picture for the weakly interacting case When uo\ ^ u>2 but uo\ rs co>2, 
a weak coupling to the COM degrees of freedom is generated, with two important 
consequences: (i) the degeneracies of the bound-state energies are lifted, and (ii) the 
coupling to the other higher-lying bound states is non-zero. Since the operators (e 
and (e commute with the ^-component L z of the angular momentum, the bound states 
are still labeled by the quantum numbers {n, m}. As far as the scattering solutions 
are concerned, the incoming wave is coupled only to states with angular momentum 
quantum number m = 0. Since (0, 0|Ce|0, 0) ~ (1, 0|Ce|1, 0), a small off-diagonal element 
(0, 0|Cs|l, 0) is sufficient to couple the bound state with {n, m} = {1, 0} to the incoming 
wave, yielding an additional CIR. As far as bound states are concerned, solutions with 
E > Eq and m = leak into the open channel, and cannot be regarded as localized 
bound states. Hence, for \a^/a\ ^> 1 and a < 0, there is only one bound state with zero 
angular momentum. In the opposite dimer limit, however, we encounter many dimer 
bound states. 

5.2. The case uo\ ^ 002 ■ Relation to experiments 

The case uu% 7^ u>2 is relevant for experiments involving two different atom species trapped 
in magnetic or optical traps El EH] For instance, in optical traps the confining 
potential depends on the detuning A = u; las — hc/X of the laser frequency co>i as from 
the characteristic frequency hc/X associated with the optical transition ns np, and 
is therefore different for two different atom species. This conclusion also applies to 
magnetic traps if the atoms are confined in hyperfine states with different projection of 
the magnetic moment along the magnetic field. As a concrete example, let us consider 




(66) 




(67) 
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a mixture of Bosonic 87 Rb atoms and Fermionic 40 K atoms. Sympathetic cooling has 
allowed to create an ultracold mixture of these two elements. By loading such a gas 
into a dipole trap and sweeping an external magnetic field, it has been possible |T3] to 
identify three heteronuclear Feshbach resonances and to measure the 3D interspecies 
scattering length a = — 14 nm. It seems feasible to tune the magnetic field near a 
Feshbach resonance and to observe the interspecies CIR. It is hence very interesting to 
know how many of them can be expected and to study their locations. 

The confining potential for a neutral atom in a standing optical wave E(r, t) = 
E (r)Re[exp (-io^J)] is Konf(r) = -(£ /4)a'E 2 (r), where a' = -e 2 / ' (2m e uj Xas E A) 
is the real part of the polarizability [Hj. Let us consider a red-detuned laser field 
corresponding to A < and a' > 0. In this configuration, the atoms are trapped around 
the maximum of the electric field. For a mixture of two species, each species experiences 
its own detuning Ak (A^b) given by the two transition wavelengths Ak = 767 nm and 
Aftb — 780 nm. Within a parabolic approximation for the potential around its minimum, 
the ratio a>K/W,b of trap frequencies for K and Rb atoms becomes 



Let us estimate this ratio for typical parameters. In order to suppress spontaneous 
emission, we assume an average detuning of A = (Ak + ARb)/2 = — 0.1co>i as , yielding 



w las = 5/ic(A K 1 + A Rb )/ll and A Rb /A k = (5A^ - 6A Rb )/(5A Rb - 6A^ ) = 0.84. Taking 



also into account the mass ratio m Rb /m K = 87/40, we have ojK/^Rb = 1-35, indicating 
a substantial coupling of COM and relative degrees of freedom. 

Using Eq. (|58|L we can project the Green's function G 4 (Rj_, 0; R' ± , 0) on the 
appropriate basis defined in Eq. (}3*Tj) and then compute numerically C,e by performing 
the imaginary-time integration, see Appendix B. Then (e can be diagonalized, and 
the effective interspecies ID interaction constant guy follows according to Eq. (|46|). 
The results are shown in the upper viewgraph of Fig. El in terms of the characteristic 
length a M = y2/G"(^K + ^Rt>))- We find two resonances, indicating that the discussion 
of Sec. 15.1.31 applies to this particular case. In order to illustrate the interpretation 
of the CIR in terms of Feshbach-type resonances with bound states of the closed 
channels, we also plot in the lower viewgraph of Fig. Elthe dimensionless binding energy 
Q = 2(E — ljk + co>Rb)/(wK + <^Rt>) of the corresponding bound state. As expected, the 
resonances occur at those values of a^/a for which the energy of the bound state of the 
closed channels coincides with the continuum threshold of the open channel. 

6. Non-parabolic confining potentials 

Describing the potential created by an optical or a magnetic guide as parabolic is to 
some extent a simplification which has to be verified. In fact, even though the lower- 
energy transverse states can rather well be approximated by the eigenfunctions of a 2D 
harmonic oscillator, in every real trap the confinement is to some degree non-parabolic. 
For resonant scattering, we expect to have a virtual occupation of many non-parabolic 




(68) 
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transverse states. As a consequence, the location of the CIR will be slightly moved, and 
new resonances could be created. This can already be seen from an analysis similar to 
the one in Sec. 15. 1.31 for small non-parabolic corrections. In order to tackle the problem 
quantitatively, a full numerical treatment is required since no analytical expression for 
the Green's function is in general available, in contrast to Sec. El 

As an example, we consider the small non-parabolicity due to the presence of 
a longitudinal magnetic bias field B z in a magnetic waveguide containing a single- 
species gas. This is necessary to avoid Majorana spin flips [THj and the subsequent 
escape of atoms out of the trap. A magnetic trapping potential is formed according to 
Konf(x) = /i m |B(x)|, where B(x) is the applied magnetic field and \i m = mpgFf'B, with 
mp being the magnetic quantum number of the atom in the hyperfine state \F, w,f), Qf 
the Lande factor and hb the Bohr magneton. Assuming that apart from the longitudinal 
bias field, the remaining magnetic fields create a parabolic and isotropic confinement in 
the transverse direction, the total confinement is given by 



where we have scaled energy in units of the parabolic trapping frequency u and length 
in units of = (/itu) -1 / 2 . The parameter x — w/ {HmB z ) is related to the Majorana spin 
flip rate T\ oss [TH] . The ID effective interaction strength #id can be calculated following 
our general approach. We compute (e numerically as outlined in Appendix B. The 
results are shown in Fig. El for x — 0.067, which corresponds to T\ oss = 10~ 6 u;. We find 
two resonances reflecting the cylindrical symmetry of the potential ()69|) and the weakness 
of non-parabolic corrections. The degeneracy of the parabolic case (shown in Fig. |3j is 
lifted and the original CIR is split into two nearby resonances. As expected, the effect of 
the non-parabolic transverse states shows up only in the deep resonant region, making 
the parabolic solution a very good approximation away from the resonant region. In 
turn, this requires a good experimental resolution in order to observe the two CIR. 

7. Conclusions 

To conclude, we have presented the general solution for two-body s-wave scattering in 
a two-component ultracold atom gas longitudinally confined to one dimension by an 
arbitrary trapping potential. The underlying key property is that the center-of-mass 
and the relative degrees of freedom of the two-particle problem do not decouple, as it is 
the case for a one-component gas and a pure parabolic confinement. Thus, no reduction 
to an effective single-particle problem is possible and the full coupled system has to 
be solved. In the framework of the pseudopotential approach, we derive the energy 
of the bound state when all transverse channels are closed. Simple analytical results 
were obtained in the limiting cases of the dimer as well as the BCS limit. Moreover, 
scattering solutions have been obtained when just one transverse channel is open. The 
effective ID interaction constant gm can be calculated after diagonalizing a reduced 
Green's function. This can be achieved analytically for the special case of parabolic 




(69) 
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confinement, where the well-known confinement-induced resonance is recovered. For 
a two-component gas, as well as for a non-parabolic confinement, more than one CIR 
occur, which reflect the symmetry properties of the confining potential. These findings 
were illustrated by applying our formalism to experimentally relevant questions. We 
are confident that once the CIR has been verified experimentally, also the effects of a 
non-parabolic trapping potential will be discerned. 
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Appendix A: Short-time Green's function 



In this Appendix, we illustrate how to expand the Green's function 

G t (X; X') = <X| exp [-(K(U) + U{X))t) |X') (70) 

with respect to t yielding the expression in Eq. (J2~3~j) for G t (Rj_, 0; R' ± , 0). In order to 
simplify the notation, we have introduced the five-dimensional vectors X = {R^,r} 
and II = {P±,p} and the functions K(U) = P\/2M + p 2 /2/i and U(X) = Vi(R± + 
/ir_ L /m 1 )-|-V2(Rj_ — /ir_|_/m 2 ) for the kinetic and the potential energy, respectively. First, 
we expand the Green's function around the free solution given by 



(X|exp [-K(U)t]\X'} 



M 
2wt 



exp 



(R ± - R' ± ) 2 M 
2t 



/ fj, \3/2 


(r — r') 2 /i" 




2t 



(71) 



o- 



In order to justify such an expansion, note that for t 
(X\K(U) exp [-K(U)t]\X') = — — <X| exp [-K(U)t]\X') oc 6(X - X" 1 



dt 



t 



(72) 
(73) 



whereas 

{X\U(X)exp [-K(U)t]\X') = U(X)(X\exp[-K(U)t}\X'} cx U(X)5(X - X') . 

Since the kinetic energy in Eq. ()72j) diverges whereas the potential energy in Eq. (|75j) 
remains finite, the latter can be regarded as a small perturbation. This expansion yields 



G t (X; X') ~ (1 - tU{X)){X\ exp [-K(U)t]\X') . 
Let us now set X = {R' ± , 0} in Eq. (|7T|) and expand with respect to t: 
(R ± , 0| exp [-K(n)t]\R! ± , 0) 



(74) 



\27TtJ 

\27TtJ 
(—) 

\27TtJ 



3/2 


M 




exp 


3/2 


fd 2 P ± 




J (2tt) 2 


3/2 


fd 2 P ± 




J (27T) 2 


3/2 






(s(R ± 



(R± - R' ± ) 2 M 
2t 

exp [iP ± ■ (R ± - R' ± 



p2 

1 - =±t ) exp [iP. 



2M 



2M 
[R± - R ± )] 



R', 



2M 
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In the last line, the operator V\ stands for (2tt)- 2 / d 2 P_ L (R L |P_ L )P3_(P_ L |R/ L ). 
Inserting X into Eq. ()74|1. we finally obtain Eq. (|23p. 

Appendix B: Evaluation of the operators Qe and £e 

In this Appendix, we outline the evaluation of the kernels (e and (e given in Eqs. (|2^|) 
and (jnOJ), respectively. 



5i; Parabolic confinement, U\ 7^ co> 2 

First, let us consider the special case of parabolic confinement, but the two species 
may experience different trap frequencies. For this confinement, the Green's function 
G t (Rj_, 0; R'j_, 0) is given in Eq. (}58j) . The first step is to project this operator onto 
the appropriate orthonormal basis defined in Eq. f!31|) . Note that this definition 

allows an arbitrary choice of the basis, apart from properly fixing the vector |0). One 
possibility is introduced in Eq. (JHEJ)- This is a natural option because it reflects the 
cylindrical symmetry of the problem. However, this choice would not permit further 
analytical progress. For this reason, we employ an alternative basis defined by 

(R±\j) = (R±_\n x ,n y ) = — t ■„ 



X 



*i>T> 



y 

CLm 



(75) 



where ip n (x) is the eigenfunction for the ID oscillator in dimensionless units, ip n (x) = 
{\fTi2 nx n\)~ 1 / 2 exp(— x 2 /2)H n (x), with H n (x) being Hermite polynomials. Note that 
the x and y directions factorize in the Green's function ()58|) . allowing to perform the 
x and y integrals separately. For convenience, we introduce dimensionless coordinates 
x — > x/a M and find 

\G t (R ± , 0; R' ± , 0)\m x , m y ) 



H 13(1-0) 



-UJlt 



with 



[F{t)]n,, 



2nt 



dxdx'ipn (x) exp 



-2ui t 



it 1 



x 2 + x' 2 



[F(t)] nx>mx [F(t)] ny>my (76) 



f(t)+xx'g(t) 



1pm (x' 



(77) 



The functions f(t) and g(t) are defined in Eq. (joTJj) . We perform the first integration by 
using the identity [TUj 



dze- {z - z ' )2 HJaz) = tt 1 / 2 (1 - 



a 2 Y' 2 H n 



az 



[l-a 2 y/\ 

with a = a(t) = [(1 + f(t))/2]~ 1 / 2 , z = x/a(t) and z' = g(t)a(t)x' /2, yielding 
[F(t)] n , m = (2 n+m m\n\)- 1 / 2 a(t)(l - a(t) 2 ) n ' 2 



X 



dx' exp 



J2 



-x ' 



a- 2 {t) 



g(t)a 2 (t) 



4 



g(t)a 2 (t) 

2(1 -« 2 (t))V2 



X 



H m (x ) . 



(79) 
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By substituting Eq. (|76|) together with Eq. (J7HJ) into (e defined in Eq. ()22|) . and by 
introducing the dimensionless time t! = y/t(uji + uty, we get 



Aixa 



x 



— / dt! \ Ah E (t') 
r<V Jo 



t 



12 



t 



12 



U)\ + Cg> 2 



2 



10) 



with the dimensionless parameter A = 2tc 3//2 /5(l — (3)a 2 Ja\ { and 



exp 


(lui + u 2 - £)f 2 " 


^1 — exp 


" 2^ " 


+ cu 2 




-( 


1 — exp 


" 2u; 2 t /2 " 


r 




UJi + cu 2 . 







(81) 

It is now possible to evaluate the matrix elements of [Cs]n,m by numerically computing 
the double integrals in Eq. (|8()J) . Note that the integrand does not suffer from any 
singularity due to the rescaling of the integration variable. Moreover, the convergence 
of the x' integral (|79|) is exponentially fast. The first term in the integrand of the t' 
integral decays exponentially at large times. Hence for large times, only the second 
term yields a contribution, where the integration can be performed analytically in this 
region. For the case of interspecies scattering of Rb and K in an optical trap, all the 
parameters entering in A,fiE{t) and [F(t)] njm can be expressed in terms of the ratios 
^Rb/^K and Aftb/AK- The generalization to determine (e is straightforward and not 
detailed further. 



B2: Non-parabolic confinement 

A numerical evaluation of the operator ( E and (e is less straightforward when the 
Green's function G t (Rj_, 0; R' ± , 0) cannot be computed analytically. In this case, 
Gi(R_i_, 0; Rj_, 0) should be computed by numerical diagonalization of the Hj_ ti and 
inserting their eigenvalues and eigenf unctions into Eq. (|2U|). For large t, this is feasible 
because only a small number of eigenfunctions contribute to the sum. However, for 
t — > 0, the number of eigenvectors required to cancel the divergence in Eq. (|22|) quickly 
proliferates. This practical limitation can fortunately be circumvented by the following 
trick. Let us formally rewrite Eq. (J22j) as 

r°° dt r 

( E (R±,R' ± )= J -{e^[G,(R ± ,0;R' ± ,0)-G t °(R ± ,0;R' ± ,0)] 

+ e Et G° t (R x , 0; R' ± , 0) - 6(R ± - R' ± )} 

r°° dt 

= / ^-e £t [G t (R ± ,0;R' ± ,0)-G°(R ± ,0;R' ± ,0)] 
Jo ^t 1 

+ C£(R±,R^), (82) 
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where G° (Rj_, 0; R'j_, 0) and ^(Rj_, Rj_) are the Green's function and the integral kernel, 
respectively, for an arbitrary reference confining potential Vo(xj_). If G°(Rj_, 0; R' ± , 0) 
is known analytically, we can deal with £g(Rj_, Rj_) as in the previous section. For 
confining potentials close to the parabolic case, we choose a parabolic Vo(xj_). 

Regarding Eq. (|82)1. we proceed as follows. We restrict the infinite-dimensional 
Hilbert space to the Af lowest eigenstates of the potential Vo(xj_), and diagonalize the 
original Hamiltonian in this AA-dimensional Hilbert space. With the eigenf unctions at 
hand, the Green's function can be computed using Eq. (|2(J|) . Then, the sum in Eq. 
(120}) is exchanged with the t-integration and the latter is performed. Next, we project 
the Green's function onto a known single-particle basis To achieve numerical 

convergence, we increase the Hilbert space dimension M until the result does not change 
anymore. We emphasize that the overall result converges to the exact result although 
obviously not all the single-particle states used in computing the Green's function 
are reliable on very long distances (comparable to the numerical system size) because 
higher-lying energy states are increasingly inaccurate. Nevertheless, the central part (in 
position space) of the eigenfunctions - which corresponds to the kinetic energy and does 
not feel the confinement - is accurate enough to cancel the divergence stemming from 
the kinetic part. In order to compute the scattering solution, we compute (e with an 
analogous procedure, diagonalize (e numerically, and insert the result into Eq. (pftjj) . 
For the non-parabolic confinement in Sec. El a parabolic Vo(x_i_) is appropriate. In this 
case, we use for {|m<)} the orthonormal basis defined in Eq. ()56|) . Then ( E is diagonal 
and given by Eq. (jHEJ). 
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Figure 1. Bound state energies E as a function oia^/a for harmonic confinement with 
equal frequency lu for both atoms. The dashed red curve indicates the bound state 
energy of the ground state |0, 0). Its binding energy Eb is given by the distance to the 
horizontal dashed line indicating the continuum threshold for the open channel. The 
blue curve marks the bound state energy of the virtual bound state relevant for the 
low-energy scattering. It is obtained by a vertical shift of the ground-state energy by 
2u, and coincides with the three- fold degenerate bound-state energy indicated as solid 
curve. The black dotted curves give the bound-state energies of the excited transverse 
states, obtained by a vertical shift of the ground-state result. 



Confinement-induced resonances in arbitrary quasi-lD traps 



22 




Figure 2. Upper viewgraph: CIR in the effective interspecies ID interaction constant 
(?id as a function of a^/a. We consider a two-component atom gas of 40 K and 87 Rb, 
with average detuning A = — 0.1o;i as (solid line). For comparison, we also show the 
result for the case when the two species experience the same trap frequency (dashed 
line). Lower viewgraph: Dimcnsionlcss binding energy f2 for the two states |0, 0) and 
|1,0>. 




Figure 3. ID effective interaction strength gin for the non-parabolic potential of Eq. 
(|69|l (black solid line) for \ — 0.067 corresponding to ri oss = 10~ 6 u;. For comparison, 
the parabolic case is also shown (red dashed line). 



